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TOWARD  A REVISED  NBS  HANDBOOK 
OF  MATHEMATICAL  FUNCTIONS 

DANIEL  W.  LOZIER 


Abstract.  A modernized  and  updated  revision  of  Abramowitz  and  Stegun, 
Handbook  of  Mathematical  Functions,  first  published  in  1964  by  the  National 
Bureau  of  Standards,  is  being  planned  for  publication  on  the  World  Wide  Web. 
The  authoritative  status  of  the  original  will  be  preserved  by  enlisting  the  aid  of 
qualified  mathematicians  and  scientists.  The  practical  emphasis  on  formulas, 
graphs  and  numerical  evaluation  will  be  extended  by  providing  an  interactive 
capability  to  permit  generation  of  tables  and  graphs  on  demand.  This  report 
provides  background  information  and  early  status  of  the  project. 


1.  Introduction 

Most  would  agree  that  the  National  Bureau  of  Standards  Handbook  of  Mathe- 
matical Functions  has  been  a phenomenal  success  in  the  field  of  mathematical  pub- 
lishing. Issued  in  June  of  1964  by  the  US  Government  Printing  Office,  this  large 
hardback  volume  of  more  than  1000  pages  has  been  in  print  continuously  since  that 
date  as  a government  publication^ . Over  200  copies  were  sold  in  a recent  one-year 
period.  A paperback  version  has  been  photo-offset  and  sold  by  Dover  Publications 
since  1965.  Although  Dover  does  not  reveal  sales  data,  it  undoubtedly  outsells  the 
government  edition  many  times  over.  Both  complete  and  abridged  editions  (with 
many  of  the  numerical  tables  removed)  have  been  produced  by  other  commercial 
publishers,  and  translations  have  appeared  in  languages  other  than  English  (for 
example,  Russian).  The  handbook  is  cited  very  frequently  in  scientific  articles  that 
make  use  of  special  functions,  making  it  one  of  the  most  cited  of  all  mathematical 
reference  works. 

What  were  the  reasons  for  this  success?  Four  might  be  comprehensiveness,  au- 
thoritativeness, timeliness  and  applicability. 

First,  comprehensiveness.  A concentrated  effort  was  made  to  identify  the  most 
important  mathematical  functions  as  measured  by  their  value  in  fields  outside 
mathematics,  such  as  physics  and  statistics,  and  to  assemble  the  most  important 
properties  of  these  functions  together  with  algorithms,  tables  of  numerical  values. 


Key  words  and  phrases.  Handbooks  of  mathematics,  special  functions,  electronic  publication. 
World  Wide  Web. 

Comments  on  any  aspect  of  a revised  Handbook  of  Mathematical  Functions  will  be  grate- 
fully received.  Further  information  and  an  e-mail  address  for  providing  comments  is  available  at 
http://math.nist.gov/DigitalMathLib. 

This  paper  was  presented  at  the  SIAM  Annual  Meeting  at  Stanford  University  on  July  14, 
1997,  in  the  Minisymposium  on  Handbooks  for  Special  Functions  and  the  World  Wide  Web.  The 
Minisymposium  was  sponsored  by  the  SIAM  Activity  Group  on  Orthogonal  Polynomials  and 
Special  Functions. 

^For  current  information  enter  the  stock  number  “003-003-00279-8”  (with  the  double-quote 
characters  included)  into  the  search  field  of  http://www.access.gpo.gov/su_docs/sale/prf/prf.html. 
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graphical  representations  and  other  relevant  information.  Most  of  the  assembled 
information  had  been  published  previously  but  it  was  unified  and  presented  in  a 
consistent  format.  A serious  effort  was  made  to  include  the  latest  results  from  the 
research  literature.  In  the  case  of  tables  of  numerical  data,  accuracy  wels  verified 
and  gaps  were  filled  by  new  computations. 

Authoritativeness  was  achieved  by  establishing  a supervisory  scientific  committee 
with  Philip  M.  Morse  of  the  Massachusetts  Institute  of  Technology  as  its  head.  The 
committee  members  were  A.  Erdelyi,  M.C.  Gray,  N.C.  Metropolis,  J.B.  Rosser,  H.C. 
Thacher,  Jr.,  John  Todd,  C.B.  Tompkins  and  J.W.  Tukey.  The  idea  to  produce 
such  a reference  work  originated  with  Milton  Abramowitz  of  NBS.  He  developed 
the  overall  plan  and  led  the  selection  of  the  authors  for  the  various  chapters.  Irene 
Stegun,  also  of  NBS,  took  over  as  technical  editor  after  the  death  of  Abramowitz. 
One  sign  of  the  authoritative  nature  of  the  handbook  is  the  impact  it  has  had 
on  standardizing  definitions  and  notation.  All  mathematical  functions  are  subject 
to  definitional  and  notational  variations  that  can  be  confusing  when  articles  by 
different  authors  are  compared,  especially  when  different  countries  or  disciplines 
are  involved.  By  becoming,  in  effect,  the  standard  reference  for  special  functions  in 
much  of  the  world,  the  handbook  helped  reduce  this  problem. 

The  appearance  of  a compendium  such  as  the  NBS  Handbook  was  timely  in  that 
it  appeared  at  just  about  the  time  electronic  computers  were  beginning  to  make 
their  influence  strongly  felt  in  applied  mathematics.  Thus  it  served  to  sum  up  the 
state  of  a highly  developed  mathematical  technology  at  a time  of  transition  to  a 
new  technology.  In  the  earlier  technology  mathematical  problems  were  considered 
solved  if  results  could  be  obtained  in  terms  of  elementary  or  higher  functions.  Be- 
cause the  properties  and  behavior  of  the  functions  were  understood,  insight  into  the 
nature  of  the  solution  was  available.  This  insight  could  be  qualitative  or  quantita- 
tive. Both  kinds  of  insight  remain  equally  valuable  today.  Qualitative  insight  comes 
from  viewing  the  ‘forest’  at  the  expense  of  the  ‘trees.’  A view  of  the  ‘forest’  coming, 
for  example,  from  knowledge  of  the  asymptotic  behavior  or  distribution  of  zeros  of 
functions  entering  into  a solution  can  be  immensely  valuable  in  scientific  and  en- 
gineering applications.  On  the  other  hand,  in  detailed  mathematical  models  of 
physical  processes  or  engineering  structures  the  ‘trees,’  e.g.  floating-point  numbers, 
are  unavoidable.  Solutions  obtained  by  more  direct  methods  such  as  quadrature, 
finite  differences  and  finite  elements  provide  an  immense  amount  of  quantitative 
information,  and  their  role  in  scientific  computing  hais  aiscended  with  the  increasing 
power  of  computers  since  1964.  However,  these  methods  provide  no  qualitative 
insight,  and  competent  practitioners  still  recognize  the  importance  of  special  func- 
tions. Three  areas  can  be  listed  in  which  they  are  important  today.  First,  they 
provide  a source  of  test  problems  for  direct  solvers.  Second,  they  are  used  to  make 
a computation  more  efficient  in  some  applications  when  length  and  time  scales  vary 
over  many  orders  of  magnitude  through  the  technique  known  as  subgrid  modeling. 
And  third,  they  can  be  highly  effective  when  solutions  are  representable  in  a series 
of  special  functions  or  as  an  integral  transform. 

The  final  reason  for  the  success  of  the  NBS  Handbook  is  its  emphasis  on  appli- 
cability. The  functions  chosen  for  inclusion  were  determined  by  their  importance  in 
applied  mathematics,  physical  sciences,  engineering  and  statistics,  both  as  a vehicle 
for  computation  and  as  a means  for  intuitive  understanding.  The  approximations 


TOWARD  A REVISED  NBS  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS 


3 


and  numerical  methods,  including  the  now  obsolete  method  of  interpolation  in  ta- 
bles, were  oriented  toward  practical  concerns,  including  the  construction  and  test- 
ing of  mathematical  software.  Accordingly,  a highly  mathematical  style  with  formal 
definitions,  theorems  and  proofs  was  avoided  while  an  insistence  on  mathematical 
rigor  was  maintained.  Lengthy  explanations  and  mathematical  developments  were 
similarly  avoided.  Thus  the  style  is  terse  and  assumes  considerable  mathematical 
capability  on  the  part  of  the  user.  The  style  and  choice  of  content  undoubtedly 
contributed  to  the  broad  and  continuing  appeal  of  the  handbook. 

The  name  we  have  adopted  for  the  envisioned  successor  to  the  NBS  Handbook  of 
Mathematical  Functions  is  the  NIST  Digital  Library  of  Mathematical  Functions^ . 
Our  approach  is  to  retain  the  ingredients  of  the  original  success  while  taking  ad- 
vantage of  a new  technological  transition.  This  transition  is  characterized  (in  part) 
by  the  use  of  powerful  computers  and  network  communications  to  disseminate 
highly  technical  reference  information  electronically.  Electronic  publishing,  online 
databases,  and  generation  of  technical  data  on  demand  are  emerging  components 
of  this  transition.  Through  its  Mathematical  and  Computational  Sciences  Division 
in  conjunction  with  its  Electron  and  Optical  Physics  Division  and  its  Standard  Ref- 
erence Data  Program,  NIST  can  provide  the  necessary  leadership  and  some  of  the 
necessary  technical,  computational  and  financial  resources.  But  outside  technical 
and  financial  assistance  is  essential,  and  toward  this  end  an  invitational  workshop 
was  held  July  28-31, 1997,  at  NIST.  The  remainder  of  this  report  takes  a closer  look 
at  the  existing  handbook  and  examines  a few  of  the  issues  involved  in  developing  a 
modernized  and  updated  version. 

2.  Existing  Scope  and  Content 

The  29  chapters  of  the  existing  handbook  are  listed  in  Table  1 . Except  for  a few 
chapters  devoted  to  introduction,  methodology  or  application,  a typical  chapter 
addresses  a particular  class  of  special  functions. 

Five  chapters  are  introductory.  Chapters  1 and  2 provide  mathematical  and 
physical  constants,  and  conversion  factors  between  metric  and  US  customary  units. 
Chapter  3 covers  binomials,  progressions,  and  means;  inequalities;  formulas  from 
calculus  and  complex  analysis;  absolute  and  relative  errors;  infinite  series;  solution 
of  quadratic,  cubic  and  quartic  equations;  successive  approximation  methods;  and 
continued  fractions.  Chapter  4 covers  the  logarithmic,  exponential,  circular,  in- 
verse circular,  hyperbolic  and  inverse  hyperbolic  functions.  Chapter  28  give  basic 
information  on  computer  arithmetic. 

Numerical  analysis  is  the  subject  of  the  methodological  Chapter  25.  It  covers  for- 
ward, central,  mean,  divided  and  reciprocal  differences;  Lagrange,  Newton,  trigono- 
metric and  inverse  interpolation,  including  formulats  with  throwback  and  stencils 
for  bivariate  interpolation;  differentiation  formulas,  including  stencils  for  partial 
derivatives  and  the  Laplacian  and  biharmonic  operators;  integration  formulas,  in- 
cluding open  and  closed  Newton-Cotes,  Gaussian,  and  other  types  of  quadrature, 
and  numerous  stencils  for  multidimensional  integration;  and  numerical  solution 
of  ordinary  differential  equations,  including  Runge-Kutta  and  predictor-corrector 
methods.  For  nearly  all  formulas  in  Chapter  25,  the  corresponding  error  order, 
error  estimate  or  error  bound  is  given. 

^In  1988  the  National  Bureau  of  Standards  became  the  National  Institute  of  Standards  and 
Technology. 
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Chapter 

Title 

1 

Mathematical  Constants 

2 

Physical  Constants  and  Conversion  Factors 

3 

Elementary  Analytical  Methods 

4 

Elementary  Transcendental  Functions 

5 

Exponential  Integral  and  Related  Functions 

6 

Gamma  Function  and  Related  Functions 

7 

Error  Function  and  Fresnel  Integrals 

8 

Legendre  Functions 

9 

Bessel  Functions  of  Integer  Order 

10 

Bessel  Functions  of  Fractional  Order 

11 

Integrals  of  Bessel  Functions 

12 

Struve  Functions  and  Related  Functions 

13 

Confluent  Hypergeometric  Functions 

14 

Coulomb  Wave  Functions 

15 

Hypergeometric  Functions 

16 

Jacobian  Elliptic  Functions  and  Theta  Functions 

17 

Elliptic  Integrals 

18 

Weierstrass  Elliptic  and  Related  Functions 

19 

Parabolic  Cylinder  Functions 

20 

Mathieu  Functions 

21 

Spheroidal  Wave  Functions 

22 

Orthogonal  Polynomials 

23 

Bernoulli  and  Euler  Polynomials,  Riemann  Zeta  Function 

24 

Combinatorial  Analysis 

25 

Numerical  Interpolation,  Differentiation  and  Integration 

26 

Probability  Functions 

27 

Miscellaneous  Functions 

28 

Scales  of  Notation 

29 

Laplace  Transforms 

Table  1.  Chapters  of  the  NBS  Handbook  of  Mathematical  Functions. 


Three  chapters  are  oriented  toward  applications.  Chapter  24  addresses  combi- 
natorial analysis  (binomial  and  multinomial  coefficients;  partitions  of  integers;  and 
Mobius  and  Euler  functions).  Chapter  26  deals  with  statistics  (definitions  and  prop- 
erties of  distribution  functions;  normal  and  bivariate  normal  probability  functions; 
and  chi-square,  incomplete  beta,  variance-ratio,  and  Student’s  t-distributions).  Fi- 
nally, chapter  29  concerns  Laplace  transforms  (definitions  and  formulas;  and  short 
tables  of  Laplace  and  Laplace-Stieltjes  transforms). 
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The  remaining  20  chapters  form  the  core  material  of  the  handbook.  Each  treats 
an  individual  class  of  special  functions,  and  typically  each  is  divided  into  four 
sections.  The  first  section,  Mathematical  Properties,  presents  definitions,  differen- 
tial equations,  integral  representations,  recurrence  relations,  functional  relations, 
series  expansions,  continued  fractions,  asymptotic  expansions,  polynomial  approx- 
imations, special  values,  derivative  formulas,  integrals,  zeros,  graphical  representa- 
tions, and  so  on.  The  second  section.  Numerical  Methods,  gives  advice  on  how  to 
compute  function  values  effectively  by  using  numerical  tables,  interpolation  meth- 
ods, and  relevant  mathematical  properties  such  as  recurrence  relations.  The  third 
section,  References,  is  divided  into  two  parts,  one  for  texts  and  the  other  for  tables. 
The  final  section,  Tables,  contains  data  taken  from  previously  published  sources, 
verified  and  augmented  with  computations  performed  at  NBS. 

3.  Opportunities  and  Obstacles 

In  the  introduction  to  this  report  an  allusion  was  made  to  an  unfolding  revo- 
lution in  information  technology  which  is  marking  a new  technological  transition. 
This  revolution  is  being  driven  by  technological  advances  in  hardware,  software,  and 
especially  in  digital  communications.  This  report  is  not  the  place  to  try  to  define 
or  describe  this  revolution  in  general  terms.  Everyone  is  well  aware  of  its  existence 
and  of  its  potential  for  dramatic  impacts  in  book  publishing,  journal  publishing,  in- 
terpersonal communication,  and  dissemination  of  all  kinds  of  information.  Indeed, 
the  impact  has  already  been  enormous.  In  mathematical  publishing,  LATEX^  is 
used  routinely  to  typeset  arbitrarily  complex  material.  Diskettes  and  compact  disks 
(CD-ROMs)  are  distributed  with  published  books  to  provide  relevant  computer  files 
to  the  purchaser.  Automated  repositories  are  used  to  obtain  software  documents 
and  other  material  by  a simple  electronic  request;  the  requester  is  the  only  human 
in  the  transaction.  These  are  examples  of  one-way  communication  from  an  author 
to  a reader  or  ‘information  consumer’  in  which  new  technology  results  in  an  im- 
proved product.  The  information  transferred  (mathematical  text,  graphs,  tables 
of  numbers,  programs)  is  prepared  in  advance  in  the  expectation  that  it  will  meet 
future  needs,  and  of  course  there  are  always  unanticipated  needs  that  are  not  met. 

What  is  needed  to  meet  unanticipated  needs  is  a mechanism  that  can  generate 
and  deliver  information  on  demand.  This  is  achievable  using  the  World  Wide  Web 
(WWW)  in  conjunction  with  computer  programs.  In  a typical  transaction,  the 
consumer  enters  a request  for  information  on  a ‘question  page,’  and  an  ‘information 
provider’  returns  an  ‘answer  page.’  The  answer  page  could  be  a simple  image  of  a 
record  extracted  from  a database  (reminiscent  of  one-way  communication  between 
an  author  and  a reader),  or  it  could  be  generated  by  an  arbitrarily  complicated 
computation.  To  get  beyond  the  traditional  concept  of  a fixed  and  immutable 
publication,  the  term  digital  library'^  has  come  into  use. 

There  are  opportunities  and  obstacles  in  applying  digital  library  technology  to 
the  subject  matter  of  the  NBS  Handbook.  The  overall  opportunity  is  to  provide 
to  the  widest  possible  audience  an  authoritative  and  up-to-date  resource  center  for 
special  functions.  We  envision  constructing  a resource  center,  or  digital  library,  in 
two  major  stages.  The  first  stage  will  be  similar  in  concept  to  the  existing  hand- 
book with  as  many  extensions  as  possible  in  the  direction  of  dynamic  generation  of 


^L.  Lamport,  LATEX  Document  Preparation  System,  Second  Ed.,  Addison-Wesley,  1994. 
■^For  example,  see  http://www.nsf.gov/od/lpa/news/publicat/frontier/5-97/5digital.htm. 
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information.  A large  database  of  static  answer  pages  will  be  constructed  and  used 
as  the  basis  for  satisfying  information  requests.  The  same  database  could  serve  as 
the  basis  for  books  and  CD-ROMs.  An  immediate  advantage  of  CD-ROMs  and  the 
WWW  over  books  is  the  possibility  of  downloading  crucial  artifacts,  such  as  sym- 
bolic code  for  mathematical  formulas,  into  documents  and  computer  programs.  In 
this  initial  stage  a major  obstacle  will  be  in  connection  with  the  numerical  tables. 
First,  the  precision  and  range  of  the  existing  tables  are  inadequate  to  current  needs, 
particularly  when  used  to  validate  numerical  software.  Second,  software  packages 
exist  that  can  compute  vastly  extended  numerical  tables — but  the  state  of  the  art 
in  software  packages  today  is  such  that  reliability  of  the  computed  values  remains 
open  to  question.  The  decision  of  what  to  include  in  the  way  of  numerical  tables,  if 
any,  will  be  difficult.  References  without  endorsement  to  algorithms  and  software 
can  be  given  and  will  be  useful  but  a completely  satisfactory  outcome  with  respect 
to  the  numerical  tables  is  not  likely  in  the  initial  stage. 

In  the  second  stage  we  envision  making  a serious  attempt  to  overcome  the  numer- 
ical tables  obstacle.  The  idea  is  to  develop  a standard  reference  library  of  numerical 
software  that  is  capable  of  computing  function  values  to  at  least  quadruple  preci- 
sion over  very  large  ranges  of  the  input  arguments  and  parameters.  But  we  require 
even  more — that  the  accuracy  of  quadruple  precision  is  reliable  in  the  sense  that  all 
errors  arising  from  truncation  and  rounding  are  monitored  and  kept  under  strict 
control.  This  requirement  is  necessary  to  uphold  the  standard  of  authoritativeness 
that  is  desirable  for  a genuine  successor  to  the  NBS  Handbook.  The  library  will  be 
used  to  generate  numerical  tables  on  demand. 

The  obstacles  to  the  construction  of  such  a library  are  mathematical  and  com- 
putational. Algorithms  with  strict  bounds  on  truncation  and  rounding  errors  are 
not  generally  available  for  special  functions.  The  techniques  of  interval  analysis  are 
applicable  here  but  most  of  the  effort  in  that  field  has  been  in  other  directions. 
The  fiaws  in  computer  arithmetic  systems,  such  as  their  vulnerability  to  underfiow 
and  overflow,  make  it  difficult  to  construct  robust  and  reliable  software.  But  these 
obstacles  provide  an  opportunity  for  creative  mathematicians  and  computer  scien- 
tists to  contribute  toward  the  solution  of  an  important  problem:  replacing  standard 
tables  with  standard  software  of  comparable  quality. 

The  use  of  standard  tables  to  validate  computer  software  was  mentioned  earlier. 
A large  effort  to  develop  advanced  test  procedures  for  mathematical  software  is 
underway  in  the  Information  Technology  Laboratory  at  NIST.  One  component 
of  this  effort  is  a Software  Test  Service  for  Special  Functions^.  This  service  will 
use  the  WWW  to  help  software  users  and  developers  craft  detailed  tests  based  on 
comparison  to  a standard.  The  most  important  missing  piece  is  a standard  reference 
library  for  special  functions.  Stopgap  meaisures  are  being  used  until  such  software 
becomes  available.  As  the  second  stage  of  our  anticipated  project  progresses,  the 
test  service  is  expected  to  become  an  integral  part  of  the  NIST  Digital  Library  of 
Mathematical  Functions. 

Graphics,  applications  of  special  functions,  and  the  role  of  symbolic  processing 
are  important  also.  For  graphics,  as  with  tables,  dynamic  generation  of  information 
promises  useful  benefits.  For  example,  fixed  graphs  could  be  replaced  by  graphical 
objects  that  could  be  manipulated  (rotated,  enlarged,  viewed  in  cross  sections,  etc.) 


®For  further  information  see  http://math.nist.gov/stssf. 
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under  the  control  of  commands  given  by  the  user.  In  the  area  of  applications,  an  ex- 
citing opportunity  exists  to  improve  communication  between  mathematicians  and 
scientists.  The  notation  and  normalization  conventions  used  in  different  fields  can 
lead  to  irritating  difficulties  in  understanding  as  well  as  actual  errors.  We  envision 
a set  of  application  digital  libraries  that  call  upon  the  Digital  Library  of  Mathemat- 
ical Functions  for  standard  data.  These  data  will  be  translated  into  a form  that  is 
immediately  meaningful  to  scientists  in  the  application  area.  Toward  this  end,  and 
as  a prototype  for  future  developments,  our  project  plan  includes  construction  of 
application  digital  libraries  of  electromagnetism  and  quantum  physics  by  the  Elec- 
tron and  Optical  Physics  Division  of  NIST.  Finally,  symbolic  processing  is  expected 
to  be  useful  also  in  helping  users  resolve  differences  in  notation  and  normalization, 
and  for  other  purposes.  User  controls,  as  in  the  case  of  manipulating  graphs,  are 
applicable  here  as  well. 


4.  Prospectus 

A reference  work  like  the  NBS  Handbook  can  be  made  even  more  valuable  by 
putting  it  on  the  WWW.  This  theme  has  been  explored,  with  attention  paid  to 
obstacles  as  well  as  opportunities. 

The  basic  content  of  the  handbook  is  found  in  the  Mathematical  Properties  sec- 
tions of  the  core  chapters.  A well-funded  project  to  revise  the  handbook  will  provide 
the  opportunity  to  update  these  sections  to  take  account  of  mathematical  devel- 
opments of  the  last  30-plus  years.  For  example,  the  chapter  on  elliptic  integrals 
predates  the  influential  work  of  Carlson  that  not  only  supplanted  the  standard 
functions  of  Legendre  with  more  satisfactory  symmetric  functions  but  also  yielded 
effective  algorithms  that  are  now  found  in  virtually  all  computer  libraries  of  spe- 
cial functions.  And  there  have  been  important  developments  in  discrete  orthogonal 
polynomials,  basic  and  generalized  hypergeometric  functions,  and  q-functions  that 
are  not  covered  at  all  in  the  existing  handbook.  We  plan  to  follow  the  proven  suc- 
cess of  Abramowitz’s  approach:  contract  with  key  mathematicians  who  can  do  a 
first-rate  job  of  writing  the  chapters  of  the  new  handbook,  under  the  supervision  of 
a carefully  selected  panel  of  senior  research  mathematicians  and  scientists  who  have 
a proven  record  of  accomplishment  with  respect  to  special  functions.  We  plan  to 
add  to  this  approach  a very  strong  validation  process,  supervised  by  the  same  panel, 
in  which  the  chapters  are  peer-reviewed  independently  of  the  authors.  Independent 
validation  was  not  a feature  of  the  original  handbook  project  but  was  used,  for  ex- 
ample, with  great  success  in  the  later  development  of  the  NAG  Library®.  Additional 
informal  reviews  and  comments  from  the  mathematical  and  scientific  community 
at  large  will  provide  still  further  validation. 

We  envision  a concentrated  effort  over  a 2-3  year  period,  after  which  the  initial 
stage  of  the  handbook  revision  project  will  be  complete.  This  period  will  begin  upon 
receipt  of  adequate  funding  authority.  The  initial  stage  will  cover  all  the  existing 
handbook,  with  extensions,  but  not  with  standard  reference  software  substituted  for 
standard  reference  tables.  The  second  stage  is  in  some  ways  already  started  through 
its  connection  to  the  Software  Test  Service  for  Special  Functions.  However,  much 
research  and  development  is  needed  before  an  acceptable  standard  library  for  the 


®B.  Ford  and  J.C.T.  Pool,  The  Evolving  NAG  Library  Service,  in  W.R.  Cowell  (Ed.),  Sources 
and  Development  of  Mathematical  Software,  Prentice-Hall,  1984,  pp.  375-397. 
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numerical  generation  of  special  functions  can  be  constructed.  Therefore  the  second 
stage  will  undoubtedly  continue  well  beyond  the  completion  of  the  initial  stage. 
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